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ABSTRACT 



Aims. We present the spectral analysis of a large set of XMM-Newton observations of EXO 0748-676, a bright dipping low-mass 
X-ray binary. In particular, we focus on the dipping phenomenon as a result of changes in the properties of the ionised gas close to 
the source. 

Methods. We analysed the high-resolution spectra collected with the reflection grating spectrometer on board XMM-Newton. We 
studied dipping and persistent spectra separately. We used the Epic data to constrain the broad-band continuum. We explored two 
simple geometrical scenarios for which we derived physical quantities of the absorbing material like the density, size, and mass. 
Results. We find that the continuum is absorbed by a neutral gas, and by both a collisionally (temperature T ~ 70 e V) and photoionised 
(ionisation parameter logif ~ 2.5) absorbers. Emission lines from OVII and OVIII are also detected. This is the first time that evidence 
of a collisionally ionised absorber has been found in a low-mass X-ray binary. The collisionally ionised absorber may be in the form 
of dense {n > 10''* cm"^) filaments, located at a distance r t 10" cm. During dips, the photoionised absorber significantly increases 
its column density (factor 2-4) while becoming less ionised. This strengthens the idea that the colder material of the accretion stream 
impinging the disc is passing on our line of sight during dips. In this scenario, we find that the distance from the neutron star to the 
impact region (~ 5 x 10'" cm) is similar to the size of the neutron star's Roche lobe. The gas observed during the persistent state 
may have a flattened geometry. Finally, we explore the possibility of the existence of material forming an initial, hotter portion of a 
circumbinary disc. 

Key words. Atomic processes - Stars: binaries - X-rays: binaries - X-rays: individuals: EXO 0748-676 



1. Introduction 

Low-mass X-ray binaries (LMXB) show different types of vari- 
ability in their Ught curve, including X-ray bursts, eclipses, and 
characteristic obscuring events called dips. A dip shows up in 
the light curve as a decrease in count-rate, which can be as 
deep as an eclipse, but with a much less clear start and end. 
, Dips are probably caused by the line of sight intercepting the 
' impact region of the accretion stream on the accretion disc 
(IWhite & Swan]dll982h . and are there fore seen in syste ms with 
inclinations in the range of 60° - 80° jFranket al.lll987l) . 
Two different models have bee n used to explain the dipping 
emission. In the first model (e.g. iParmar et alJll986) two com- 
ponents (both power laws with an exponential cut off) are used 
to model the dipping emission. The first component is absorbed 
and the second one is not. This approach is sometimes referred 
to as the "absorbed + unabsorbed approach". 
The second model, which consists of two types of emitting com- 
ponents, was proposed bv IChurch et al.l (Il998l) . The first com- 
ponent is point-like blackbody emission from the neutron star 
itself, and the second component is extended Comptonised emis- 
sion from the accretion disc corona. During dipping intervals, the 
accretion disc corona is progressively covered, hence the name 
"progressive covering approach". 



Recently. lBorrin et al.l (l2005l) have proposed a different approach 
to modelling the dipping spectra in the LMXB 4U 1323-62. 
They find that the spectral changes during dipping emission can 
be modelled by varying the properties of a photoionised absorber 
(i.e. column density and ionisation parameter). No absorbed 
+ unabsorbed components or progressive and partial covering 
of an extended corona is necessary and abs orption lines and 
the co ntinuum are modelled self consistently. IDfaz Trigo et alj 
( l2006h show that this approach can be used on a number of other 
X-ray binaries exhibiting dipping behaviour, among which also 
EXO 0748-676. 

EXO 0748-676 is an LMXB that shows all the variabil- 
ity mentioned ab ove: burst s, eclips es, and dips. It was dis- 
covered in 1986 ( Parmar et a l. 1986: Gottwald et al.lll986l) and 



has been stu died regularly s ince its disc overy ( Parm ar et al 
1991: Hertz et alJll995L[T997l:lThomas et al,, 1 997: Church et al 



19981: ISidohetal.1 120051: IWolff etalj l2005h especially since 
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the launch of the high resolution X-ray observatories XMM- 
Newton (Bonnet-Bidaud et al. 2001: Cottametal. 2001, 2001 
Homan et al. 2003: Diaz Trigo et al.l 120061) . and Chandra 
(Jimenez-Garate et al..,2003i) . 

iParmar et alj (Il986l) derived the basic properties for EXO 
0748-676. The LMXB has a period of 3.82 hours as measured 
by X-ray eclipse timings, and the eclipses themselves have a du- 
ration of 8.3 minut es. Assuming a compact object of 1.4Mo, 
IParmar et al.l (Il986h infer the mass of the companion to be 
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between O.OSSMq and O.45M0 and the system inclination be- 
tween 75° and 82°. Recently a strong X-ray burst was observed 
(Wolff et al. 2005) from which a distance to EXO 0748-676 was 
derived that depends on whether the burst is hydrogen dominated 
(5.9 + 0.9 kpc) or helium dominated (7.7 + 0.9 kpc). In this paper 
we use the average value of 6.8 kpc. 

We study here the ionisation state of the accretion disc 
around the neutron star of EXO 0748-676 in order to get a bet- 
ter understanding of the geometry, dynamics and ion stratifica- 
tion of the disc. We use all RGS data available in the wavelength 
range 7-35 A. 

In Section |2] we describe the observations and the data anal- 
ysis, in Section[3]we present the results of the spectral fitting. In 
section |4] we discuss the results and derive physical parameters 
of the system and in Section|5]we present our conclusions. 

2. Data analysis 

The data were taken from the 'KMM-Newton arch ives and were 
acquired with the Epic PN (IStriider et alj 1200 Ih. Epi c MOS 
dTurner et al.' 2001) and RGS (Iden Herder et all BOOll ) instru- 
ments. For the analysis of each data set we require first the avail- 
ability of the RGS and then of either the PN or MOS cameras 
(Table[Tll- We processed the data with version 6.5.0 of the XMM- 
Newton Science Analysis System (SAS). For the Epic analy- 
sis, the background was extracted from a circular region in the 
same field of view of the source. From this region we produce 
a light curve from which we removed any background flares. 
Only four observations were totally discarded because of a per- 
sistently flaring background. To obtain the source light curve and 
spectrum we use an annular extraction region (7?,„ 8" and 
Rout ~ 51.5") to prevent pile-up. We divide the source radiation 
into hard (5 < E(keV) < 10) and soft (0.3 < E(keV) < 5) X-ray 
emission. We use the hard X-ray light curve to obtain the starting 
and ending times for the bursts and the e chpses, since this energ y 
range is marginally affected by dipping jPiaz Trigo et al.ll2006h . 
EXO 0748- 676 occasion ally shows two, or even three, bursts in 
a row ( Boirin et al.ll200 7'). with an inter-burst time that is much 
shorter than between two normal bursts. We count these double 
and triple bursts as one, since we are only interested in the start- 
ing and ending times of the total interval containing bursts. 
We also use the PN data to obtain the Good Time Intervals (GTI) 
for the dipping and eclipsing emission (we will call persistent 
emission everything that is not dip, burst or eclipse). These GTIs 
were subsequently used to extract the RGS spectra. If, for a cer- 
tain observation, PN (or MOS) were not switched on at the same 
time as RGS, we only use the overlapping intervals. 
Next, we divide the hard X-ray light curve by the soft X-ray light 
curve. We use the resulting hardness ratio curve to distinguish 
between dipping and persistent emission: a dip interval shows 
up as an increase in the hardness ratio. We use the average of 
a representative portion of the hardness ratio curve that clearly 
belongs to the persistent emission to select dipping and persis- 
tent emission: values higher than 2 times the average persistent 
hardness ratio were selected as dipping emission. 
Upon close inspection of the hardness ratio curves we found that 
it was not possible to distinguish between persistent and dipping 
emission in the way that was described above for the observa- 
tions from orbits 0692 and 0708. The reason for this is that the 
hardness ratio curve for these observations does not show clear 
separate periods of dipping and persistent emission as compared 
to a normal case. For this reason, even if these two data sets sat- 
isfied our first selection criteria, we did not analyse them any 
further. 



We process the RGS data with the SAS task "rgsproc". After we 
check the source coordinates and correct them if necessary, we 
create a spectrum of the RGS. 

We apply an effective area correction to the RGS data processed 
with SAS version 6.5.0 (private communication of J. Vinkfl We 
fit the area-co rrected spectra wit h the X-ray and UV fitting pack- 
age SPE?fl ( Kaastra et alj|1996l) . The best fit was found by min- 
imisation of tiiex^- Errors are at Icr unless otherwise stated. 



2.1. The data set 

In table [T] we list the observations from the XMM-Newton 
archive used in this paper. In the first column the orbit in which 
the observation was made is listed. The second column shows 
the observation identifier and the third column shows which in- 
struments were used to observe the source. The "obs. date" col- 
umn gives the date on which the observation was made. The la- 
bel "observing time" gives the time subdivided in "total", per- 
sistent ("per") and dipping ("dip") times. Differences are caused 
mainly by the fact that eclipses and bursts are left out. 
We grouped the observations, when possible, (column "group" 
in table [TJ using the SAS-task "rgscombine". The first criterion 
for combining separate spectra is to look whether the observa- 
tions were close in time. Secondly, we fit a continuum model 
including absorption (as defined in the next sections) to the data 
and calculate the errors of the model parameters. If the parame- 
ters of the separate observations are consistent within the errors, 
we combine the observations. 

We combine the observations from orbits 0055, 0067 and 0212 
into group A, the observations from orbit 0338 into group B and 
the observations from orbits 06XX/07XX in three groups. In par- 
ticular, we kept the observation from orbit 0710 as a separate 
group CI. We then combined the observations from orbits 0693 
and 0719 into group C2, those from orbits 0694 and 0695 into 
group C3 (Table[B. 

We use a single observation from each of the groups defined 
in Table [1] to create a spectrum of the PN (or MOS, if no PN 
is present) data and we fit both the RGS spectra and the PN 
spectrum of a group simultaneously to obtain the best fit. The 
observation we use to obtain the PN spectra is indicated in the 
last column of table [T] There are two reasons to use PN data. 
The first is to constrain the power-law photon index in the wave- 
length range of the RGS instrument. The second is that PN data 
may offer us further constrains on a possible highly ionised ab- 
sorber, whose feat ures would show up in the 6-7 keV band (e.g. 
iBoirin et alj|2005b . 



2.2. Spectral modelling 

Two representative RGS spectra of persistent and dipping emis- 
sion are plotted in Fig. [1] In both spectra we see a number of 
interesting absorption features. At 23.1 A we see an absorp- 
tion edge that is produced by neutral oxygen in the interstellar 
medium. We also see two other absorption edges, at 14.2 A and 
16.8 A due to OVIII and OVII respectively. We can also clearly 
discern two emission lines, OVII at 21.8 A and OVIII at 19 A, 
the latter being more evident in the dipping spectrum. 



more 
website of 



information can be 
tiie XMM-Newton 



found 
users 



tlie 
group: 



http://xmm.vilspa.esa.es/external/xmm_user_support/usersgroup/20060518/rgs_calib. 
Tliis correction lias been incorporated into subsequent SAS releases. 
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orbit 


obs. ID 


instruments 


obs. date 


observing time (ks) 


group 


spectrum 








(yymmdd) 


total 


per 


dip 




mos/pn 


0055 


0122310301 


mos/rgs 


00-03-28 


61.0 


11.9 


0.9 


A 


- 


0067 


0123500101 


mos/pn/rgs 


00-04-21 


62.1 


16.9 


14.5 


A 


yes 


0212 


0134561101 


mos/pn/rgs 


01-02-03 


8.4 


0.0 


3.8 


A 


- 


0212 


0134561201 


mos/pn/rgs 


01-02-04 


6.6 


0.3 


1.7 


A 


- 


0212 


0134561301 


mos/pn/rgs 


01-02-04 


6.9 


1.6 


1.3 


A 


- 


0212 


0134561401 


mos/pn/rgs 


01-02-04 


6.6 


1.5 


1.4 


A 


- 


0212 


0134561501 


mos/pn/rgs 


01-02-04 


6.6 


1.2 


1.3 


A 


- 


0338 


0134562101 


mos/rgs 


01-10-13 


8.7 


4.7 


2.1 


B 


- 


0338 


0134562201 


mos/rgs 


01-10-13 


6.9 


4.8 


1.3 


B 


- 


0338 


0134562301 


mos/rgs 


01-10-13 


6.9 


6.2 


0.0 


B 




0338 


0134562401 


mos/rgs 


01-10-13 


6.9 


2.8 


2.3 


B 




0338 


0134562501 


mos/rgs 


01-10-13/14 


6.9 


2.5 


2.4 


B 


yes 


0693 


0160760201 


mos/pn/rgs 


03-09-21/22 


93.4 


41.5 


23.1 


C2 




0694 


0160760301 


mos/pn/rgs 


03-09-23/24 


108.3 


47.8 


29.6 


C3 


yes 


0695 


0160760401 


mos/pn/rgs 


03-09-25/26 


82.0 


29.4 


26.6 


C3 




0710 


0160760801 


mos/pn/rgs 


03-10-25/26 


67.9 


21.8 


18.1 


CI 


yes 


0719 


0160761301 


mos/pn/rgs 


03-11-12/13 


94.7 


31.3 


35.3 


C2 


yes 




10 15 20 25 50 35 10 15 20 25 30 35 

Wavelength (1) Wavelength (1) 



Fig. 1. Two representative spectra of RGSl and RGS2, taken from obsid 0160760201. The left plot shows the persistent spectrum 
and the right plot shows the dipping spectrum. In both plots the upper panel shows the data (points) and model (solid line) in Counts 
s ' A ' . The dotted line at the bottom is the background. The lower panel shows the model residuals ((data-model)/error). The data 
have been rebinned for clarity. 



As a starting point, we modelled the RGS spectra following 
iDiaz Trigo et al. 12006), who analysed the lower resolution data 
of Epic-PN (orbit 0719). In particular their absorption model 
consisted of the combination of a neutral absorber and one pho- 
toionised absorber The continuum was modelled by a emission 
from a black body and a powerlaw. In our case, applying the 
same model to the RGS wavelength range of 7-35 A does not 
yield a satisfactory fit. In particular, the model is not capable of 
fitting the OVll and OVlll edges simultaneously, strongly sug- 
gesting the presence of a second, ionised, absorbing component. 

We therefore departed from this starting model to fit the 
high-quality, high-resolution RGS data. For clarity, in the fol- 
lowing the modelling is referred to the data set shown in Fig.[I] 
(left panel). In addition to the neutral absorber and one ionised 
absorber (with log^ ~ 2.5, where ^ is defined as the ratio 
between the ionising luminosity and the product between the 
gas density and its squared distance from the central source: 
we included first an additional photoionised ab- 
sorber (model xabs in SPEX). This did not yield a satisfactory 
fit (x^ldof = 2276.31/1188 = 1.92, for the persistent emis- 



sion data. The column density of this additional component is 
~ 1.8 X 10^' cm"^ with ionisation parameter log^ ~ 0.8. The 
OVll edge region, close to the iron unresolved transition array 
(UTA), between 14 and 20 A, is where the model fails the most 
in fitting the data (Fig.|2]i. Next we tried a more complex absorp- 
tion model, which mimics a continuous distribution of column 
densities as a function of the ionisation parameter (model warm 
in SPEX). The warm model did not improve the fit significantly 
with respect to a single-ionisation-parameter gas component. 
The best fit is reached when a collisionally ionised gas compo- 
nent is added to the model (model hot in SPEX) with a;^'^/dof = 
1631/1188 - 1.37 for this specific data set. The gas temperature 
is ~ 75 eV and the gas column density is ~ 1.5 x 10^' cm"^. In 
Fig.|2]we compare the two interpretations in the 14-20 A region 
where the OVll edge and the iron UTA creates a broad trough 
in the data. In the upper panel the transmission due to the pho- 
toionised (light line) and the collisionally ionised gas (dark line) 
are shown. In the next panel the data and the models, convolved 
by the instrument resolution, are shown (with the same colour 
labelling as above). It is evident that the photoionised absorber 
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model causes a too shallow absorption at 16.6-17 A. At the same 
time this model severely underestimates the data at 17.5-18 A, 
where absorption by OVI and OVII takes place. On the contrary, 
the coUisionally ionised absorption model provides a good fit 
to the data. The reason is that given a temperature, a collision- 
ally ionised gas ionic column density distribution strongly peaks 
around a single ion (OVII in this case), while for a photoionised 
gas, a range of ions around a main ion is allowed. This pro- 
duces the unwanted OV and OVI absorption in the xabs model 
in Fig.|2] 

The parameters for the absorption components included in 
the final fit are listed in Table |2] The underlying emission is 
produced by a power law, a blackbody and two Gaussian to 
model the OVII triplet and OVIII Lya emission lines (Table |2]l. 
The Gaussian profiles FWHM ranges from 0.25 to 0.5 A for 
OVII and from 0.05 to 0.35 A for OVIII, respectively. The main 
component of the OVII triplet Une arises around the wavelength 
of the intercombination line at 2 1.8 A. There is no evidence for 
the forbidden line at 22.1 A. The second emission line in the fit 
is the OVIII Lya line, at 18.96 A. The value of the reduced 
reported in Table |2] reflects the cross calibration uncertainties 
between RGS and PN that can be up to 10-15%. However we 
are interested in narrow features located in a spectral range 
where instruments are well calibrated. We verified that a slight 
change of the continuum shape does not aff'ect our results. 



3. Results 

The persistent and the dipping state of EXO 0748-676 can be 
interpreted using the same model for all data groups (Table |2]i. 
While the central source flux changed over time, the continuum 
shape remained stable. On the contrary, the absorbers showed 
significant variability. In Fig. [3] we show the values of the col- 
umn density of the three absorbers and the ionisation parameter 
of the ionised absorber as a function of time. The column den- 
sity of all three absorbers systematically increases during dip- 
ping emission. 

Fig. [3] (first panel) shows that during dips, the column den- 
sity of the neutral gas increases. As we cannot disentangle any 
intrinsic neutral absorber from the interstellar one, we see this 
as an increase of the total neutral absorption. However, the col- 
umn density in the direction of the system is ~ 1.01 x 10^' cm"^ 
dKalberla et al.l l2005h . comparable to the neutral absorber in 
the persistent emission. This suggests that neutral (or mildly 
ionised) material must exist in the vicinity of the source, in par- 
ticular near the gas producing the dips. 

The coUisionally ionised absorber that models mainly the 
OVII edge (Sect. |272l Fig.|2]i generally shows an increase in col- 
umn density during dipping, while the temperature of this com- 
ponent does not show a clear change. The different behaviour 
of this ionised absorber in group A (see Table |2] and Fig. 13) is 
probably caused by the low statistics of that group in the dipping 
spectrum. 

The photoionised absorber shows the same behaviour tha t 
was already found for this source by iDfaz Trigo et al.l (l2006l) : 
during a dip the column density increases, while the ionisation 
parameter decreases. Compared to the persistent intervals the ab- 
sorbing material b ecomes denser and less ionised during dips 
jBoirin et al.ll2005l) . 

For most groups the values of the blackbody normalisation 
and temperature overlap within errors when we compare dip- 
ping and persistent emission. A moderate scattering is expected 




14 



19 



20 



15 16 17 18 

Wavelegth (A) 

Fig. 2. Comparison between the collisionally-ionised gas model 
(hot model, dark line) and a photoionised gas {xabs model, light 
line) in the OVII edge-iron UTA wavelength region. This is a 
zoom of Fig.[T] left panel. From top to bottom, first panel: the- 
oretical transmission curve. Second panel: data points together 
with the hot and xabs models. The xabs model fails to fit the 
data due to extra absorption by OV-OVI in the 17.5-18 A region 
and in the iron UTA region. Third and fourth panels: residuals in 
terms of cr to the xabs and hot models, respectively. 



due to slow changes of the persistent emission during long ob- 
servations. For the power-law parameters the values of F and 
normalisation, appear instead more scattered when we compare 
dipping and persistent emission. This may be partially due to 
cross-calibration eff'ects between RGS and PN. 



4. Discussion 

4.1. Evidence for a coUisionally ionised gas 

We find that the spectrum of EXO 0748-676 can be best mod- 
elled by adding a coUisionally ionised absorber This model well 
fits both the persistent and the dipping emission in all data sets. 
This is the first time that a coUisionally ionised absorber is found 
in a low mass X-ray binary system. The column density of all the 
absorbers present in our model increases during dips suggesting 
an association with the gas producing the dip. A coUisionally 
ionised plasma can exist in the presence of a high ionising ra- 
diation field {L ~ 2 X 10^^ erg s"') only if it is located at a suf- 
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group A group B group CI group C2 group C 3 



parameter 


pers 




dip 




pers 






dip 




pers 


dip 




pers 


dip 




pers 




dip 






neutral abs. 










































CI 


1.28 ± 


0.05 


3.81 ± 


0.8 


1.10 + 





03 


2.21 


±0.17 


1.11 ± 0.05 


2.36 


±0.13 


0.82 ± 0.02 


1.48 


±0.03 


0.84 


±0.20 


1.25 


± 


03 


4.05 ± 


0.25 


1.8 ± 1 


.6 


0.95 ± 





10 


5.28 


±0.56 


T oo+O 18 
■^■^'^-0:40 


5.9 ± 


0.6 


1.95 ± 0.07 


4.1 ±0.1 


1.09 


±0.08 


3.42 


± 


13 




66 ±3 




< 80 




61 ± 4 






63 ± 


5 


70 ± 3 


70 ± 


7 


73 ± 2 


69 ± 


1 


80 ± 


4 


65 ± 


2 




< 
PI 

"h 


< 12.27 


0. fixed 


< 4.83 






123 + 30 


< 7.6 


< 11 


1 


< 4.7 


< 1 1 


.1 


< 3.57 


< 15 


.91 




41 ±4 




160+'° 




37 ± 4 






74 ± 


11 


46 ± 4 


140 ± 10 


39 ± 4 


80 ± 


3 


19 ± 


2 


51 ± 


2 






2.35 ± 


0.04 


2.24 ± 


0.04 


2.56 ± 





04 


2.36 


±0.05 


2.40 ± 0.04 


2.39 


±0.03 


2.63 ± 0.04 


2.38 


±0.01 


2.51 


±0.04 


2.33 


±0 


02 




< 8.7 




< 6.1 




< 32.5 






21 ± 


9 


< 48 


11+^ 




17 ± 2 


13. ± 


3 


14 ± 


2 


11 ± 


2 




powsr law 












































1.07 ± 


0.05 


1.60 ± 


0.16 


l.c5U ± 


U 


Uo 


1.11 


±0.17 


i. / u ± U.U / 


1.92 


±0.12 


i .0 J ± U.Uj 


1.15 


±0.02 


1.91 


±0.04 


1.49 


±0 




p 


1.23 ± 


0.02 


1.42 ± 


0.05 







0^ 

\jj 


1.34 


±0.08 




1.46 


±0.03 


1.44 4- 0.01 


1.24 


±0.02 


1.47 


±0.01 


1.35 


±0 


01 


blackbody 










































norm'^ 


3.4 ±0.5 


22+''' 

^ -20 




1.20 + 





13 


3.6 ± 


1.0 


1.8 ±0.4 


14. + 


6 


1.10±0.12 


1.95 


±0.25 


1.10 


±0.09 


1.88 


±0 


21 




0.10 ± 


0.01 


0.08 + 


0.02 


0.14 + 





01 


0.12 


±0.01 


0.11 ±0.01 


0.08 


±0.01 


0.12 ±0.02 


0.11 


±0.01 


0.12 


±0.02 


0.11 


±0 


02 


O VII 










































norm'' 


11 ± 1 




136+*^ 




3.3 +C 


.9 




27+'^ 




9± 1 


40 ± 


8 


3.6 ±0.4 


9.6 ± 


0.8 


2.9 J 


:0.4 


7.6 ±0.7 




FWHM^ 


0.28 ± 


0.04 


0.35 ± 


0.06 


< 0.5 






22+''-2* 


0.23 ± 0.03 


0.44 


±0.07 


0.28 ± 0.05 


0.36 


±0.05 


0.24 


±0.06 


<0.1 






O VIII 










































norm*" 


2.7 ± 0.3 


16!^ 




1 6+'-'^ 






c +40 
J._2 




9 iq+3.74 
^■-•^-047 


S 40+67.96 


0.51 ±0.16 


2.32 


±0.26 


0.87 


±0.21 


1.56 


±0 


24 


FWHMf 


0.16 ± 


0.04 


0.34 ± 


0.09 


0.07 + 





07 


< 0.11 


0.07 ± 0.05 


0.05 


±0.05 


20+° '° 


10+"°^ 


0.33 


±0.09 


0.19 


±0 


05 


x'/dof 


1.29 




1.41 




1.42 






1.52 




1.38 


1.50 




1.68 


1.83 




1.59 




1.66 







Notes: 

" Column densities (A'h) are given in units of 10^' cm"^ for the neutral, collisionally ionised (CI) and photoionised (PI) absorbers. 
* The width Ci. of the absorption lines is expressed in km s"' . 
Temperatures, T, are in eV for the collisionally ionised gas and in keV for the blackbody. 

The normalisation (norm) of the powerlaw is given in 10'*'' phs"' keV"' at 1 keV. The normalisation of the blackbody is in units of 10'^ cm^. 
The normalisation of the emission lines (OVII and OVIII) are expressed in units of lO''^ ph s"' . 
^ The FWHM of the emission lines is in km s ' . 



ficiently large distance from the central source. If we put the 
gas at such a large distance (r ~ few x 10" cm) and we require 
^ = L/nr^ < 1, we find n > fewx 10'"* cm"^. From the measured 
column density we then derive a very narrow thickness of the gas 
I = Nh/h < 1.5 X 10^ cm. The geometry of this gas could then 
be a uniform shell of gas, with a thickness of less than 150 km. 
More realistically, the medium could be inhomogeneous, with a 
low filling factor, as in a smoke-like medium. 
In parallel, a further test of the physical process producing the 
absorption comes from the associated three OVII emission lines. 
These c an be in principle used as a density/temperatu re diagnos- 
tic (e.g. lPorquet & Dubaull2000l: ICottam et al.ll2001l) . The inter- 
combination (x+y) line is evident in the data while due to ab- 
sorption some flux from the w line may be suppressed. There is 
no evidence of the forbidden (z) OVII l ine at 22.1 A, sugge sting 
either a high-density photoionised gas jCottam et al.ll200l]) or a 
collisionally ionised gas. For a photoionised gas we would also 
expect radiative recombination continuum of OVII, which is not 
detected in the data. 

However, both the collisionally ionised and the log^ ~ 2.5 pho- 
toionised plasma included in our model (Table|2]i do have a 
significant amount of OVII and OVIII and they can both con- 
tribute to the emission spectrum, making difficult to disentangle 
between the two ionising process contributions. The measured 
width of the lines is however large (Sect. 12. 2] ). This suggests that 
a considerable amount of matter should be located well in th e 
accretion disc, assuming Keplerian motion ( Cottam et al.ll200i1) . 
The FWHM of the detected oxygen lines is roughly 1000 km s"' . 



Assuming a I.4M0 central object then the distance must be of 
the order of 10'° cm. Therefore the emission spectrum could be 
dominated by the photoionised gas producing broad lines, hiding 
the narrow emission lines produced by the collisionally ionised 
gas. Indeed, the inferred line width of a collisionally ionised gas 
rotating at large (r ~ few x 10" cm) distance is of the order of 
100 km s ' . The absence of the OVII z line is however consistent 
with our requirement of n > few x lO'"* cm""* for the collision- 
ally ionised gas. For the photoionised gas, the suppression of the 
z line could be due either to a high gas density or to a strong UV 
radiation field illuminating the gas. 



4.2. The geometry of the system 

Here we test two simple geometries for the absorbing/emitting 
photoionised gas. We consider the absorption dur ing dips as 
caused by photoionised gas (see iBoirin et alj|2005l for full de- 
tails). In both case 1 and 2, described below, we assume that the 
dipping gas has constant density rid and it is located at a distance 
Rd from the central source in a protuberance (the impact region) 
of the disc. 

As dips are merely caused by absorption, we further assume 
that the intrinsic lumin osity of the source rema ins constant dur- 
ing all the observations lDfaz Trigo et akl (l2006h and lBoirin et al.l 
(l2007l) . 

Fig.|4]depicts the geometries described in detail below. For both 
panels, the right hand sides describes the observer's view during 
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Fig. 3. Parameter values as a function of observation date. Diamonds refer to persistent and triangles to dipping state. From top to 
bottom, left to right: neutral absorber column density, coUisionally ionised absorber column density, photoionised absorber column 
density and ionisation parameter log ^. Error bars are at Icr level. 



the dipping state. The left hand sides shows the view during the 
persistent emission. 

Case 1 : the absorption during persistent intervals is due to a 
thin shell of gas of width x and ionisation parameter (Figure 
|4]and the caption). The shell is located at a distance Rp\. We as- 
sume X to be approximately the same as the size of the compact 
impact region that produces the dips (Fig.|4] upper-right panel). 
The rationale behind this is that the gas that produces the per- 
sistent absorption could well be gas that trails from the impact 
region into the rest of the disc. With the assumption that the in- 
trinsic luminosity does not depend on the dipping state of the 
source, we can compute; 



(1) 



where ^ is the ionisation parameter defined before, n the hydro- 
gen density, R the distance to the ionising source and the sub- 
scripts p and d refer to persistent and dipping emission respec- 
tively. 

Since we assume that the density has no radial dependence, 
we can rewrite the above for case 1 as: 



■ NH,d 



'N, 



(2) 



H.p 



Here A^h = is the hydrogen column density derived from our 
fits for the dipping {Nu.d) and persistent (A^H.p) emission, (see 



gland tableH. 

Case 2: the persistent gas is in a spherical cloud encompass- 
ing the whole neutron-star-disc system (lower plot in Figure|4]i. 
This cloud of photoionised gas with ionisation parameter has 
a radius Rp2 a density rip. The dipping gas is again a compact 
region of width x (Fig.H) lower-right panel). In this second case 
we f urther make the assumption that the disc is geometrically 
thin (iFrank et al.ll2002t iLewin et al.lll995l) . Therefore, we con- 
sider the height of the disc H ~ Q.lRdi and also H ~ x (Fig.|4]i. 
The persistent gas is sphere-shaped, therefore Nh,p = Rpi^p- For 
the dipping gas A^e.d = xnd = Hnd = Q.lRd2nd, because of the 
simple relation between H and Rd2 defined above. We can now 
substitute the values of A^h in equation ([T]i and obtain the relation 
between Rp and Rd for case 2: 



Rp2 ^ 10 -^d- Nii^d 

Rd2 ^p ■ Nh,p 



(3) 



For both case 1 and case 2, we have calculated the value of 
Rp/Rd and listed them in Table [3] The values for A^h and ^ are 
taken from Tab. |2] 

We see that in case 1 the radius of the shell is comparable 
to the distance of the impact region to the central source (Rd ~ 
Rp). This is consistent with the idea that the shell could be a 
continuation of the impact region. On the contrary, in case 2 the 
cloud of gas in the persistent state is an order of magnitude larger 
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DIPPING STATE 



S I -disk 
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Rp2 











1 disk ■■ 













Rd2 



Fig. 4. Two possible geometries: case 1 in the upper plot the gas 
producing the persistent absorption has a shell-like geometry of 
thickness x and distance Rd\ from the central source. In case two 
(lower plot) the persistent gas is a cloud-like geometry of radius 
Rp2 encompassing the whole system. The labels "NS" and "disc" 
indicate the position of the NS and the accretion disc, respec- 
tively. On the right-hand side of the drawings, on both panels, 
is the impact region of size x and distance from the central 
source. The drawn line of sights indicates which environments 
the observer intercept during the persistent and dipping state. 

Table 3. The ratio of the persistent over dipping radius Rp/Rd 
for the two geometries and for each data group (A-C3). 



A B 


CI 


C2 


C3 


easel 1.7 + 0.4 1.1 ±0.1 


1.7 ±0.2 


1.1 ±0.1 


1.3 ±0.1 


case 2 30+ 13 13 ± 3 


30 ±6 


11 ±2 


17 ± 3 



(Rp » Rd) than the distance of the impact region to the central 
source. 

4.3. A flattened geometry for the persistent gas 

For both case 1 and 2 described in the previous section, it is 
also possible to put additional constraints on the geometry of the 
persistent absorbing gas. This can be achieved us ing the OVII 
and O VIII emission lines and the absorption edges (ICottam et al.l 
12001'). As justified in Sect. 14. l[ we assume here that the emission 
lines are dominated by photoionisation. In this calculation we 
only use the OVII line and edge for the data groups where the 
associated errors on those values were the smallest. Therefore 
we excluded group B. 

We start with the equation for the line luminosity of a given 
ion i dOedahl 19991): 



dVaRR(T)i],f,A [ph/s]. 



(4) 



where n^. is the electron density, ajm is the radiative recom- 
bination coefficient for the transition, rji is the fraction of re- 



combinations that produce the line, is the fractional abun- 
dance of the ion in ionisation state i+l and A is the absolute 
abundance of the element. The radiative recombination coef- 
ficient ffss at 10 eV, which corresponds to the plasma tem- 
perature for the best-fit ionisation parameter, is aRR(ovii) - 
1.6 X 10 "cm^s-i dVerner & Ferlandl Il996h . The fraction of 
recombinations tha t produce the OVII line is rjnovii) - 0.75 
(ICottam et al.l 2001.). The absolute abundance for oxygen is 
taken from lAnders & Grevessg ([19891) : A = 8.5114 x 10 The 
emission measure is defined as: 



EM 



-'line 



aRR{T)rjifA 



(5) 



The integral in equation (|5]i can be estimated for both the shell- 
like (case 1) and the sphere-like (case 2) geometry for the per- 
sistent gas. 

For case 1 , for a shell extending from r io r + x (x <^ r) w& have 



EMshell 

nx — —. 

In fact we consider r 
sphere of gas: 

3 EAl,;pj^gfg 



(6) 



= Rp. For case 2 we simply have for a 



4n nr^ 

The value of nr^ can be derived from the ionisation parameter: 



(7) 



(8) 



Inserting equation (O, the column densities of equations ^ and 
(|7]i are now a function of the ionised absorber, the line luminos- 
ity, the emission measure and the photoionisation parameter. All 
these parameters can be evaluated from the spectral analysis. 

In Table|4]we compare the expected column density (for case 
1, nx and case 2, nr) of the gas, as derived by the line emission, 
with the one we measure (from the absorbed spectrum) for each 
data group. 

Assuming at first that the material subtends a solid angle (as seen 
from the neutron star) of An sr, the emission of the line will be 
maximised and, as a consequence, also the derived column den- 
sity of the photoionised emitting gas. If we attribute any differ- 
ence to geometrical effects, the ratio between observed/expected 
column density provides us a flattening factor (Table |4]i. From 
the Table we see that the derived column density is systemati- 
cally lower than the measured one, suggesting that the material 
does not fully extended in the vertical direction, but rather it is 
flattened. 

Incase 1 the flattening factor/' is ~ 4.5-9. In case 2 the discrep- 
ancy between observed and expected column density is some- 
what lower, ~ 1 .5 - 3 . For case 1 , the flattening facto r is in agree- 
ment with that found previously for this source ( ICottam et al.l 
2001). Converting the flattening factor into an opening angle, 
the ionised persistent gas subtends an average angle of 18° for 
case 1 and a wider angle of 54° for case 2. 

4.4. Density, size, and mass oftlie pliotoionised gas 

Using the information we gathered above, we calculated the val- 
ues for the gas density n and distance Rp (Table |4]i. For both 
geometries we use the identity: n - (nr)^/nr^. 



8 



J. van Peet et al: XMM-Newton observation of EXO 0748-676 



types of binary sys t ems (e.g. cataclysmic v ariables, novae etc., 
iDeufel et al.' 1999; Solheim & Sionj 11994 . including neutron 
stars (Muno & Mauerhan 2006). The existence of such material, 
leading possibly to the formation of a circumbinary disc 
jDubus et alj|2002l) . would have important implication s for the 
evolutions of the binary system dSpruit & Taamll2b0lb . In this 
scenario, the detected collisionally ionised gas (Sect. 14. Il l could 
also be a further extension, too far away to be photoionised, 
of the persistent emission. However, the measured temperature 
of the collisionally ionised gas is ~70 eV, which translates in a 
thermal velocity v,/,p„„ of about 140 km s"'. If the gas is located 
at a distance i 10'^ cm from the central source v,/,p„„ becomes 
comparable to the escape velocity v^sc ~ 190 km s"', making the 
gas a short-lived component, unless continuously replenished by 
the accretion disc activity. However, if the distance is ~ 10" cm, 
the escape velocity akeady becomes four times larger than 
Vihenm guaranteeing the stability of the collisionally ionised 
gas. This would favour a scenario similar to case 1, where the 
distance of the persistent and dipping gasses range between 
0.5 and 1.4x10" cm. With this geometry, the collisionally 
ionised gas could well be an extension of the persistent gas, 
perhaps constituting a first, hotter, portion of a circumbinary 
disc, which may eventually join a colder phase (T ~600K) 
which extends up to four times the distance of the accretion disc 
dMuno & Mauerhan''2006). 

In the hypothesis that the collisionally ionised gas is at a distance 
comparable to the persistent and dipping material, it is likely 
that there is interaction with the system where the accretion 
stream impacts with the disc, together with neutral material. A 
higher concentration (implying a higher column density) of the 
collisionally ionised and neutral gas in correspondence wit h the 
dipping material is then possible (e.g. iBoirin et alj I2005L and 
Sect.©. 

The average gas density is of the order of /ewx lO'^ cm"^ for 
a shell-like geometry (case 1, Table |4| and about 100 times less 
for the spherical geometry (case 2, Table |4]i. However, the com- 
plete suppression of the forbidden line in the persistent spectrum 
requires a gas density > lO'^^cm"-', unless the gas is irradiated 
by a strong ultraviolet field. Therefore, if the shape of the OVII 
triplet is dominated by density effects, we have to require that n 
is not uniform in the gas, having for example a radial gradient. 
Finally, we calculate the mass of the persistent gas enclosed both 
in a flattened shell and in a sphere using M - nnipV where M is 
the mass in the gas, n is the density, nip is the proton mass and V 
is the volume of the gas (Table|4|. 

The average distance from the neutron star to the dipping 
material Rd is simply obtained using the Rp/Rd ratio (Table [3]). 

The distance from the neutron star to the inner Lagrangian point 5. Conclusions 
(LI) varies between 6.4 - 9.7 x 10'° cm d Frank et alj|2002 '). For 

both the geometries under study, the location of the dipping In this paper we present the spectral analysis of a large XMM- 

material, where the flow of matter from the donor impinges Newton data set of EXO 0748-676 mainly focusing on the RGS 

into the accretion disc, is consistent with the distance of LI. high-energy resolution data in order to investigate in depth the 

The gas seen during the persistent phase is at approximately the gaseous environment of the source. 

same distance of the dipping point in case 1, while in case 2 the For both the persistent and dipping spectra, the best fit model is 
sphere of persistent gas should extend much further. The average a combination of a neutral absorber, a collisionally ionised ab- 
thermal velocity of the photoionised gas {v,herm ~ 77kms"') sorber, a photoionised absorber, a power law, a black body and 
is more than a factor two lower than the escape velocity at a two Gaussian emission lines at the energies of OVII intercom- 
distance of 10'^ cm (vesc ~ 190km s '). Therefore in principle bination line and the OVIII Lya. The scatter in the parameter 
the photoionised gas can survive also beyond LI. values for the different data groups gives an idea of the level of 

uncertainty of our conclusions. 

Observational evidence of cold gas extending beyond the or- Here we summarise the main physical impUcations of this anal- 
bits of the components has been found in infra-red for different ysis: 



Table 4. For each data group, we list the relevant observed and 
calculated values to obtain the flattening of the gas (Sect. |4.3l l. 
For case 1 and case 2, we also derive the persistent gas column 
density, flattening, density, distance from the central source Rp 
and its mass. We also calculate the distance of the dipping mate- 
rial Rd- 







A 


CI 


C2 


C3 


observed values 




10"" ergs/s 


1.38 


1.05 


1.06 


0.97 


log^ 


erg • cm 


2.35 


2.40 


2.63 


2.51 


ATxabs 


10^2 


4.12 


4.63 


3.97 


1.94 


fz 




0.317 


0.274 


0.132 


0.200 


Lline 


lO'*- ph/s 


11.2 


9.04 


3.63 


2.97 


calculated values 


EM 


lO^'ph/cm' 


3.46 


3.23 


2.69 


1.45 


2 

nr 


lO'-* cm-' 


6.16 


4.18 


2.48 


3.00 


case 1 












nx 


10^- cm-2 


0.45 


0.62 


0.86 


0.39 






9.23 


7.47 


4.62 


4.97 


n 


10" cm-3 


27.6 


50.4 


64.0 


12.2 


Rp 


10" cm 


1.4 


0.9 


0.6 


1.5 


Rd 


10" cm 


0.8 


0.5 


0.5 


1.1 


M 


10-'2 Mo 


1.1 


0.5 


0.2 


0.6 


case 2 












nr 


10^- cm-2 


1.34 


1.85 


2.59 


1.16 






3.08 


2.50 


1.53 


1.67 


n 


10" cm-3 


0.2 


0.5 


0.6 


0.1 


Rp 


10" cm 


14.9 


9.0 


6.2 


15.5 


Rd 


10" cm 


0.5 


0.3 


0.5 


0.9 


M 


10-'- Mo 


22.5 


8.9 


3.1 


10.9 



For case 1 we adopt x ~ 0.1 x r as the thickness of the shell. We 
obtain: 



2 



.^^.O.l^.^:l^^(M . (9) 

nxf j'Lline \ ^ I 



For case 2 we have: 



nr^ \47:aRRT]f,Aj \Lc, 



con! 
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- The gas during the dipping phase has a higher column 
density and it is less ionised than the gas in the cloud that 
permanently surrounds the system. 

This is in line with the i dea of ionised absorbers causing the 
dippi ng phenomenon (iBoirin et alj 120051; iDiaz Trigo et alj 
120061) . In this scenario the compact and colder impact 
region, where the accretion stream hits the disc, is in our 
line of sight during dipping intervals. 

- We used the information obtained from the RGS spectral 
analysis to investigate the geometry of the persistent/dipping 
gas and their distance from the central source. 

First, in the hypothesis of the persistent gas in a shell form, 
we obtain that the thickness and the radius (~ fewx 10' ' cm) 
of the shell are consistent with those of the dipping gas. The 
intuitive picture is that the persistent gas is a trailing tail of 
the dipping gas. This axial symmetric geometry should be 
flattened with an opening angle of r oughly 18°, consiste ntly 
with previous studies of this source dCottam et alj|200lb . 

In a second scenario, where a spherical cloud of persistent 
gas was considered, we find a smaller flattening, as the angle 
subtended by the gas is -54°. The persistent gas extends 
about ten times further than the distance of the dipping gas. 
The radius of the gas sphere is about 0.9 - 1.5 x 10'^ cm, 
i.e. larger than the first Lagrangian point. Photoionised 
material at such a distance might be the first portion of a 
wider circumbinary dis c, detected in infra red in at least two 
neutron stars systems jMuno & Mauerhanll2006h . 

For both geometries the gas density is not sufficiently high 
to justify the suppression of the OVII forbidden fine in 
favou r of the OVII intercombination line dPorquet & Dubaul 
120001) . 

However, assuming that the emission line is broadened 
by Keplerian motion around the central source, we infer a 
lower limit for the location of the emitting gas of ~ lO''* cm, 
implying a higher gas density (~ lO'^^cm"^). This value 
is certainly sufficient to produc e the observed OVII triplet 
ratios (Porquet & Dubaul l2000h . This suggests that there 
may be a radial dependence of the gas density. Another 
possibility is that the forbidden line in the photoionised gas 
is suppressed by a strong UV field. 

- It is the first time that a collisionally ionised absorber is 
found in a low mass X-ray binary system. In order for 
this component to survive to the strong radiation field, it 
should be located far away from the central source (r ~ 
few X 10" cm). This is supported by the evidence that both 
the temperature and the column density of such gas change 
very little during the dipping or persistent phase and seem 
to be unperturbed on very long (years) time scale. We es- 
timated that this collisionally ionised absorber should be in 
the form of dense filaments, with a very low filling factor 
In the hypothesis of a circumbinary disc, this gas might also 
be an ulterior continuation of the persistent gas. In this case 
the collisionally ionised gas should not be located beyond 
few X 10" cm, in order not to easily escape the system. This 
would favour a more compact-sized geometry for the whole 
gas environment of the neutron star, as the one described in 
case 1, where the persistent gas is a trailing tail of the dipping 
bulge, at a distance of about 10" cm. 
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